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Due to the smaller axial tilt of the Moon, interiors of 
some of the polar craters on the lunar surface never 
get sunlight and are considered to be permanently 
shadowed regions. Several of such regions are known 
to contain water-ice deposits. These regions are 
expected to show elevated values of circular polariza¬ 
tion ratio (CPR). Hence, the interiors of craters con¬ 
taining water-ice deposits are characterized by elevated 
CPR values as observed by the A-band synthetic aper¬ 
ture radar (Mini-SAR) on-board Chandrayaan-1 mis¬ 
sion of ISRO. However, elevated CPR values were also 
observed from the interiors of some non-polar craters 
and also from young, fresh polar craters. Thus, ele¬ 
vated CPR values are not a unique signature of water- 
ice deposits. Therefore, additional information related 
to geological setting and roughness patterns should 
also be considered while identifying the regions con¬ 
taining water-ice deposits. For identifying a unique 
signature of water-ice deposits, analysis of radar scat¬ 
tering mechanism in elevated CPR regions was car¬ 
ried out. Areas of elevated CPR due to double-bounce 
and surface scattering conditions were segmented and 
polarimetric, backscattering properties of diffuse scat¬ 
tered were analysed. Based on the signatures of dif¬ 
fuse scatterers and radar backscattering coefficient, a 
scattering mechanism-based algorithm was developed, 
which has the advantage in classifying regions show¬ 
ing elevated CPR due to surface and double-bounce 
scattering effects. The algorithm was then tested using 
Mini-SAR data and it was also found to be useful in 
eliminating regions of elevated CPR in fresh craters 
observed due to the double-bounce effect. 


Keywords: Circular polarization ratio, craters, lunar 
surface, water-ice deposits. 

Recent investigations of lunar polar regions using 
remote sensing data suggested the existence of water-ice 
in permanently dark regions of polar craters 1-4 . Because 
of the small tilt of the Moon’s axis with respect to the 
Sun, craters and other depressions near the poles never 
receive sunlight and remain in permanently shadowed 
zones. This condition creates cold trap zones with tempe- 
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ratures as low as 50-70°K, in which water-ice may have 
existed for at least two billion years 5 . Cometary debris 
and meteorites containing water-bearing minerals 
bombard the Moon’s surface and the water is mostly lost 
in space. It is believed that the water molecules which are 
trapped in such cold traps, can never find a way to come 
out and significant quantities of water-ice would accumu¬ 
late in geological timescale. 

Radar is an optimal instrument for detecting metre- 
thick deposits of water-ice due to its distinctive polariza¬ 
tion signature compared to the surrounding silicate 
regolith 4 . One of the methods to detect such ice deposits 
is the use of radar polarimetry which has the potential to 
provide additional information about the surface and sub¬ 
surface physical properties, and it has been often demon¬ 
strated in remote sensing studies of our solar system 6 . 
The key radar parameter for detection of water-ice is the 
circular polarization ratio (CPR), which is the ratio of the 
same-sense circular (SC) echoes to the opposite-sense- 
circular (OC) echoes. In case of circular transmit polari¬ 
zation and linear received polarization system, it is the ra¬ 
tio of left circular polarized to opposite polarized signal, 
that is right circular polarization. Specular echoes from 
surfaces that are smooth at wavelength scales will lead 
to low ratios and hence low CPR values, while diffuse 
scattering from rough surface generates CPR values 
approaching one or even greater than one for extremely 
rough terrain (caused by multiple double-bounce reflec¬ 
tions). Also, low temperature water-ice when it is present 
as ice is transparent to RF energy and radar is multiple 
scattered by imperfections and inclusions in the ice 1 . Fur¬ 
ther, as observed from some planetary bodies, the ice 
inhomogeneities exhibit a phenomenon known as coher¬ 
ent backscatter opposition effect 7 (CBOE), which causes 
an increase in the radar echo reflectivity and enhance¬ 
ment of CPR along the backscatter direction. 

The Clementine bi-static radar experiment detected 
elevated CPR values in the south pole of the Moon, sug¬ 
gesting the presence of patchy ice deposits in the perma¬ 
nently shadowed areas of Shackleton crater 8,9 . This was 
later questioned by further studies using ground-based 
Arecibo observatory and Green Bank Telescope (GBT) at 
a bi-static angle of 0.37°, by suggesting that high CPR 
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values observed near the south pole are caused by surface 
roughness of the rugged inner walls of Shackleton crater 
and proximal ejecta of impact craters 10 . However, with 
the availability of data from the first orbital Synthetic 
Aperture Radar (SAR) around the Moon, the Mini-SAR 
on-board Chandrayaan-1, Spudis et al. 1 showed the 
anomalous behaviour of craters which are under perma¬ 
nent Sun shadow, i.e. elevated CPR values in their interi¬ 
ors, but not exterior to their rims. Spudis et al . 1 analysed 
CPR distributions at both the interior and exterior of a 
small polar anomalous crater, a polar fresh impact crater 
and another small non-polar anomalous crater. It was 
found that the CPR distributions of non-polar anomalous 
crater, which is exposed to full sunlight over the course 
of a lunar day, and the fresh impact crater are similar, but 
not the polar anomalous crater. Statistical analysis pre¬ 
dicted deposits of water-ice inside the polar anomalous 
crater, which correlate with the proposed locations of 
water-ice modelled on the basis of Lunar Prospector epi¬ 
thermal neutron data 11 . Thomson et al 4 suggested an up¬ 
per limit of -5-10 wt% H 2 0 ice (up to 30 vol%) present 
in the uppermost metre of regolith within the Shackleton 
crater floor at the lunar south pole. Modelling of radar 
scattering at two different wavelengths by Thompson et 
al n indicated that scattering from near-surface ice cov¬ 
ered by a thin regolith layer can be separated from rocks 
if the SC enhancement is twice the average or more. It 
was also observed that if the lunar ice is dispersed 
throughout the regolith as ice-filling pores, then scatter¬ 
ing differences might be too small to detect 12 . 

The high CPR values, thought to be due to scattering 
from water-ice, were also found to be within fresh, young 
craters due to roughness/double-bounce effect associated 
with them. The anomalous scattering properties were also 
observed at many non-polar craters on the Moon. Thus, 
CPR values greater than one alone could not be attributed 
to ice deposits. Hence for the unambiguous detection of 
water-ice, it is required to separate the high CPR areas 
caused by double-bounce effect, and the exact locations 
within the interior of the polar anomalous craters which 
contribute to volume/diffuse scattering must be identified 
and mapped. In order to eliminate regions of high CPR 
due to double-bounce, an algorithm was proposed, which 
removes the contribution of multiple double-bounce re¬ 
flections from the total signal received from areas having 
CPR values greater than one. In this case, since the SC 
signals are more enhanced than the OC signals received, 
surface scattering (Bragg scattering) contribution is neg¬ 
ligible. This algorithm was tested using Chandrayaan-1 
Mini-SAR data. 

Mini-SAR data and derived polarimetric 
parameters 

Mini-SAR imaging radar on-board Chandrayaan-1 of 
ISRO was the first lunar orbiting mono-static SAR with a 


unique hybrid polarity architecture 13 , flown in October 
2008 with a primary scientific objective to detect water- 
ice in the permanently shadowed regions on the lunar 
poles up to a depth of a few metres. Mini-SAR operated 
at the S-band (2.38 GHz) frequency with left circular 
polarization (LCP) transmission and reception in linear 
horizontal (//) and vertical ( V) polarizations. During its 
operational period, Mini-SAR mapped over 95% of the 
lunar poles (lat. >80°) at 35° incidence angle with a 
ground resolution of 150 m. This radar used Stokes’ 
parameters and derived daughter products such as CPR, 
degree of linear polarization and relative phase to 
describe the backscattered field, in order to distinguish 
volume scattering from other scattering mechanisms 14 
(e.g. sub wavelength-scale surface roughness). 

A typical Mini-SAR image strip contains 16 bytes of 
data in four channels of four bytes per pixel, given by 
\LH\ 2 , \LV\ 2 , Real (LH LV *) and imaginary (LHLV*). The 
first two channels represent the intensity images for hori¬ 
zontal and vertical receive and the last two represent real 
and imaginary components of the complex cross-power 
intensity between the horizontal and vertical receive, 
respectively. From these data, Stokes’ vector can be gen¬ 
erated 13 to completely describe the polarization state of 
the received wave as: 

~S 0 =(\LH\ 2 +\LV\ 2 ) 

S,=(\LH\ 2 -\LV\ 2 ) 

S= 1 ' 1 ' ' , (1) 

5 2 = 29T (LH-LV*) 

5 3 = -23 {LH-LV*) 

where the terms LH and LV are complex numbers denot¬ 
ing magnitudes of electric fields corresponding to LCP 
transmission and H and V receive respectively; () repre¬ 
sents ensemble averaging in time or spatial domain, * 
represents the conjugate and 91 and 3 are the real and 
imaginary parts, respectively, of the complex number. 
The sign on S 3 is negative, which is consistent with the 
backscattering alignment (BSA) convention. The Mini- 
SAR data were found to have a LH-LV relative phase ( S) 
shift of 45° based on the analysis of different portions of 
the lunar surface 15 . In order to overcome this systematic 
error, a phase calibration was done and phase offset of 
45° was applied to the Mini-SAR data before further 
processing, as explained in recent studies 15,16 . 

Several useful quantitative measures follow from 
Stokes vector 13 ; some of the important derived measures 
have been explained in Mohan et al. 15 . The decomposi¬ 
tion of m-delta (m - degree of polarization, delta - 
relative phase) feature space offers a novel way of deter¬ 
mining the odd, even bounce and volume scattering con¬ 
tributions corresponding to a target or within a radar 
resolution cell 13 . The degree of polarization has long been 
recognized as the single most important parameter char¬ 
acteristic of a partially polarized EM field, which 
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Table 1 . Average values of total backscattered power, circular polarization ratio (CPR), relative LH and LV phase (delta) and Degree of polarization 
(DoP) at interior (int) and exterior (ext) regions of various craters studied here. The CPR values indicate mean +(-) standard deviation 


Total backscattered power (dB) CPR Delta (degree) DoP (m) 


Crater diameter (km) 

Int 

Ext 

Int 

Ext 

Int 

Ext 

Int 

Ext 

Polar anomalous craters 

16 

-13.32 

-13.42 

0.7±0.27 

0.3410.13 

-60.21 

-72.31 

0.4 

0.64 

6 

-11.36 

-10.86 

0.77 10.28 

0.44 1 0.2 

-49.7 

-72.25 

0.36 

0.56 

8 

-12.7 

-13.48 

0.85 10.32 

0.37 10.15 

-38.3 

-81.85 

0.31 

0.58 

20 

-14 

-10.35 

0.7 10.26 

0.4610.16 

-65.08 

-75.61 

0.38 

0.53 

11 

-14 

-13.44 

0.57 10.21 

0.32 10.13 

-72.41 

-75.92 

0.44 

0.65 

5 

-13.53 

-12.9 

0.68 10.25 

0.3410.12 

-74.22 

-85.78 

0.37 

0.61 

4 

-13 

-12.82 

0.68 10.27 

0.43 1 0.22 

-76.82 

-81.44 

0.39 

0.53 

4 

-9 

-12.37 

0.77 10.24 

0.2410.11 

-52.74 

-84.53 

0.36 

0.74 

4 

-13 

-14 

0.74 10.31 

0.12 10.06 

-69.65 

-97.83 

0.35 

0.83 

3 

-10.86 

-9.1 

0.63 10.23 

0.22 10.08 

-77.51 

-79.83 

0.41 

0.79 

3 

-15 

-12.74 

0.55 10.19 

0.2010.06 

-95.32 

-98.4 

0.41 

0.73 

14 

-11.8 

-9.45 

0.6610.25 

0.3 10.18 

-72.73 

-83.29 

0.39 

0.58 

8 

-13.33 

-14.3 

0.84 10.31 

0.31 10.12 

-43.2 

-78.76 

0.32 

0.66 

Polar fresh craters 

20 

-12.33 

-9.26 

0.77 10.27 

0.7 10.23 

-56.2 

-59.91 

0.35 

0.39 

14 

-9.6 

-5.34 

0.81 10.29 

0.84 10.27 

-43.88 

-40 

0.36 

0.36 

3 

-13 

-3.65 

0.75 10.24 

0.71 10.26 

-57.53 

-58.81 

0.32 

0.4 

Non-polar anomalous craters 

5 

-10.31 

-9.32 

1.11 10.35 

0.610.19 

13.34 

-80.15 

0.28 

0.42 

8 

-9.16 

-11 

1.08 10.39 

0.65 1 0.26 

4.29 

-53.28 

0.31 

0.42 

3 

-11.13 

-11.32 

0.68 10.21 

0.4 10.13 

-60.6 

-71.33 

0.4 

0.6 

1 

-9.95 

-9.04 

0.85 10.29 

0.4610.14 

-31.74 

-84.44 

0.38 

0.52 



Figure 1 . North polar anomalous craters: Total backscattered power (So) and circular polarization ratio (CPR) images for each crater are repre¬ 
sented by (i) and (ii) respectively, (a) Mouchez M crater, (b) Hermite A crater (c) Small crater on the floor of the bigger crater Hermite, (d) 
Whipple crater to the north of Peary crater, indicated by an arrow. Scale bar at the bottom indicates the range of CPR values. 
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Figure 2. Non-polar anomalous craters: The So and CPR images for each crater are represented by (i) and 
(ii) respectively, (a) Unnamed crater to the north of Byrgius crater indicated by an arrow, ( b ) Kunowsky C crater, (c) 
Small crater in the Mare Orientale basin towards the south of crater Maunder, indicated by an arrow. ( d) Unnamed cra¬ 
ter on the ring of Pyroclastic Orientale deposits. Scale bar on the bottom right indicates the range of CPR values. 


represents the fully polarized case (m = 1) and the fully 
depolarized case (m = 0), and is not affected by polarimetric 
basis transformations 13 . The LH-LV relative phase ( 8) is 
highly sensitive to polarimetric variations in the backscat- 
tered field, and under the condition of circularly polarized 
illumination, it is the prime indicator of double-bounce 
scattering. A decomposition technique based on the para¬ 
meters m and 8 derived from the hybrid polarimetric SAR 
data was found to be effective in identifying dominant scat¬ 
tering mechanisms associated with a radar resolution cell 13 . 
The m-8 scatter plot gives overall information of the scat¬ 
tering contributions within the area under consideration; so 
a decomposition colour code 17 was introduced given by 


R - x m x (1 + sin 8)12 

indicates double-bounce backscatter 
(dihedral reflectors), 

G = yJS 0 x (1 -m) indicates volume scattering V 
(randomly polarized targets), 

B = \[Sq xmx (1-sin 8)12 

indicates single bounce backscatter 
(specular reflection), J 


( 2 ) 


where the total backscattered power S$ = R 2 + B 2 + G 2 . 


A circular wave backscattered by a surface rotates the 
incident wave in the opposite direction. Consequently, for 
the case of left circularly polarized transmission, a domi¬ 
nant surface scattering mechanism will return a right 
circular polarized wave (delta = -90°). So in terms of 
m-delta representation, the surface scattering term con¬ 
sists of (1 - sin£) term. This colour code was applied for 
the study regions to identify the scattering mechanism. 

Study area and methodology 

In this study, scattering properties of 20 craters were ana¬ 
lysed; 16 of them are at the north pole and the remaining 
craters belong to non-polar regions of the Moon. Out of 
the 16 north polar craters, 13 have permanently shadowed 
regions 1 (which makes up to one-third of the total north 
pole anomalous craters), while the remaining three are 
fresh impact craters with their ejecta fields distinctly 
visible in the Mini-SAR images. An inventory of craters 
analysed in this study is summarized in Table 1. The total 
backscattered power and CPR images of four polar 
anomalous craters are shown in Figure 1. To begin with, 
the CPR distributions over the interior and exterior por¬ 
tions of the north polar craters (both fresh and anoma¬ 
lous) were analysed and found to be in good correlation 
with earlier results 1 . 
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Figure 3. m-Delta decomposition images of the north polar craters: (a) Small anomalous crater on the floor of 
Rozhdestvensky. A typical double-bounce feature is indicated with an arrow. ( b ) Fresh impact crater Main L. 
(c) Whipple crater towards top of the image. Small secondary craters on the floor of Peary crater can also be seen 
towards south of Whipple. Decomposed image represents double-bounce (red), volume (green) and surface (blue) 
scattering contributions. 


A few non-polar craters were identified that exhibit 
anomalous scattering properties, as observed by Spudis 
et al.\ like the crater to the north of Byrgius (24°S, 
65.3°W; 5 km diameter). The backscatter properties of 
this crater along with three other non-polar anomalous 
craters with elevated CPRs in their interiors (Figure 2) 
were analysed. These craters are in the sunlit region 
where water-ice cannot exist and are listed from Table 1. 
The Mini-SAR data were analysed in terms of backscatter 
and derived Stokes parameters. Radar backscatter was 
analysed for the effect of surface roughness, whereas 
polarimetric parameters were used for identifying the 
scattering mechanism. 

Analysis of scattering mechanism in high CPR 
regions 

The high CPR values can be associated with areas of 
water-ice deposits, or due to very rough surfaces such as 
a rough, blocky lava flow, which has angles that form 
many small comer reflectors. In this case, the radar signal 
could hit a rock face (changing LCP into RCP) and then 
bounce over to another rock face (changing RCP back 
into LCP) and hence to the receiver, which produces a 
double-bounce effect that could mimic the enhanced CPR 
obtained from ice targets 18 . Therefore, elevated CPR 
alone cannot be considered as a derived indicator of the 
presence of water-ice particles distributed in the regolith, 
and we require to understand the scattering mechanism 
within the elevated CPR regions as well. 

Analysis of decomposition images: From the colour- 
coded decomposition images, various scattering mecha- 
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nisms present in a given region could be identified. Fig¬ 
ure 3 shows the m-delta decomposition images of three 
polar craters. The 8 km diameter crater on the floor of 
Rozhdestvensky (Figure 3 a) and the Whipple crater 
(Figure 3 c) are anomalous and the Main L crater (14 km 
diameter, Figure 3 b) is a fresh crater. As observed from 
Figure 3, the dominant scattering mechanism is volume 
scattering in the interior of both the anomalous craters 
and a mixture of even and odd bounce scattering beyond 
their rims. A characteristic double-bounce feature, as 
explained by Raney et al. 19 could be observed towards 
the western edge of the small anomalous crater on the 
floor of Rozhdestvensky (Figure 3 a). In the case of fresh 
impact crater Main L, a mixture of double-bounce and 
volume scattering (Figure 3 b) was observed to be domi¬ 
nant both inside and outside the rim of the feature. This is 
consistent with its high CPR values due to enhanced 
degree of wavelength-scale roughness caused by the 
abundant blocks present within the interior and at the 
ejecta deposits of the crater. Each of the craters along 
with its impact ejecta region was analysed in terms of ra¬ 
dar backscatter and scattering mechanism to study the 
anomalous scattering properties. Table 1 shows the aver¬ 
age values of various polarimetric parameters associated 
with the analysed craters which are So, CPR, m and 8. 
The craters are classified into anomalous and fresh 
(young) according to the mean CPR values observed. 
Since high CPR regions were to be analysed, the areas 
with CPR values greater than 1 at the crater interior were 
masked out for subsequent analysis. 

Analysis of scatter plots: The scatter plots of scattering 
mechanisms derived from m-8 decomposition method 

1583 








RESEARCH ARTICLES 



0 01 02 03 0-4 (j oi oa 0 3 0.4 03 06 0 7 09 0 01 02 03 04 0 5 0 6 

Scattering Contribution! Scattering Contribution Scattering Contribution 


Figure 4. Scatter plots of CPR versus scattering mechanism at the interior of (a) small polar anomalous crater on the floor of Rozhdestvensky, ( b) 
polar fresh impact crater Main L and (c) non-polar anomalous crater towards the north of Byrgius crater. 
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Figure 5. CPR vs delta scatter plots for the interior of three of the studied craters: (a) Small crater on the floor of Rozhdestvensky. ( b ) Fresh im¬ 
pact crater Main L. ( c ) Whipple crater. CPR is on the X-axis and delta is on the 7-axis. The colour spectrum represents the density of the pixels 
with red representing the maximum. 


CPR-deita segmentation 



Circular polarization ralio 

Figure 6. Graphical illustration of CPR-delta segmentation tech¬ 
nique. Zones of CPR less than and greater than 1 are indicated. The red 
dotted lines indicate regions of S= -180° to 0° and 5= 90° to 100° in 
those zones respectively, which are segregated using this method. The 
residual regions contribute exclusively to diffuse/volume scattering in 
the CPR > 1 zone. 
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against CPR were analysed for each class of crater as 
described in Table 1. The high CPR regions (CPR> 1) 
showed contributions from surface, volume and double¬ 
bounce scattering mechanisms for all the craters. This 
could be observed in the plots of scattering mechanism 
versus CPR (Figure 4) for the three craters shown in Fig¬ 
ure 3. Thus, high CPR regions exclusively contributed 
by diffuse/volume scattering were identified using an 
algorithm proposed here. 

CPR-delta-backscatter algorithm for water-ice 
detection 

Since the relative phase delta ( S) is a sensitive indicator 
for odd and even bounce scattering, regions showing sur¬ 
face scattering (specular or mirror-like reflections) could 
be separated from other regions. Therefore, use of CPR in 
conjunction with 8 could give improved understanding of 
scattering types in a region. For this analysis, CPR values 
were plotted against 8 (Figure 5) and the scatter plots 
were studied. In CPR space, the image pixels for which 
CPR < 1 were found to be confined to the region of 8= 
-180° to 8= 0°, and the upper half of the plane ( 8= 0° to 
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Figure 7. The CPR-delta algorithm for identification of water-ice in lunar polar regions. 



Figure 8. Regions of volume scattering (in green) overlaid on So images of the small crater on the floor 
of Rozhdestvensky using (a) existing method (regions of CPR >1) and ( b ) the algorithm developed. 
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Figure 9. Regions of volume scattering (in green) overlaid on S 0 images of Whipple crater region using 
(a) existing method (regions of CPR >1) and (b) the algorithm developed. 


8= 180°) contained pixels with CPR > 1, which have 
contributions from both double-bounce and volume scat¬ 
tering. In this way, contribution from surface scattering 
was eliminated. As the Mini-SAR transmitted LCP, a 
dominant double-bounce scattering mechanism will 
return a left circular polarized wave again, which will be 
centred at delta = + 90°. So the pixels with delta values of 
90° ±10° were removed to separate the regions having 
contribution from both double-bounce and volume scat¬ 
tering. A graphical representation of this segmentation 
using CPR and delta values is shown in Figure 6. 

Also, as a threshold value of - 15 dB was chosen based 
on the observations of mean backscatter at regions of 
CPR > 1 at the interior of the analysed craters. This back¬ 
scatter threshold was applied to mask out two regions: 
interior of non-polar anomalous craters which mimic the 
water-ice-like condition and very high backscatter 
observed at the crater walls which are oriented in the 
direction of the radar (in the order of + 5 dB to -8 dB). 
Based on the above inputs, an algorithm was developed as 
shown in Figure 7, which can be used to represent the 
areas of possible locations of water-ice. 

Results and discussion 

Table 1 gives an estimate of the average scattering pro¬ 
perties at the interior and exterior of the craters analysed 
in this study. It can be observed from Table 1 that the 
interior and exterior regions of the polar anomalous cra- 
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ters are characterized by volume scattering with low m 
values and surface scattering with relatively high m val¬ 
ues respectively. Regions both interior and exterior to the 
rim of young, fresh craters were observed to have low m 
values. In case of non-polar anomalous craters, the values 
indicate the presence of mixed scattering contributions 
with lower m values. 

The results of the existing decomposition technique 
and the algorithm developed are shown in Figures 8 and 
9. The green pixels overlaid on S 0 images indicate the 
possible locations of water-ice according to the two 
methods. Clear differences could be observed between 
the two methods in depicting the possible water-ice 
regions within the permanently shadowed regions of the 
crater interior. Even though some non-polar craters also 
exhibit anomalous scattering properties, they are charac¬ 
terized by high backscatter values in their interiors (aver¬ 
age values around -10 dB) in contrast to the low values 
of polar anomalous craters, suggestive of diffuse scatter¬ 
ing likely due to water-ice particles mixed with the 
regolith. Thus, CPR-delta method was implemented for 
detecting regions of water-ice along with a backscatter 
criterion that helps in differentiating polar anomalous 
craters from those of non-polar origin. 

Conclusion 

The CPR-delta approach along with backscatter analysis 
used in this study provides a better method to identify 
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volume scattering regions and the algorithm developed 
showed promising results for identifying craters with 
water-ice deposits by separating them from young, fresh 
craters with elevated CPR values. Identification of regions 
which exclusively contribute to volume scattering in the 
permanently shadowed regions is important in supporting 
the debate over the possible existence of water-ice in 
those areas, since a high degree of correlation between 
the occurrence of high CPR values and the transition 
from seasonally sunlit to permanently shadowed regions 
of the crater wall is certainly required. The algorithm 
showed promising results for identifying craters with 
water-ice deposits and separating them from young, fresh 
craters with elevated CPR values. 
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